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I. INTRODUCTION 



Nonequilibrium quantum field theory became a topic of interest as a result of several cru- 
cial questions facing physics in the 1980's. Before that time, people were mainly concerned 
with calculating cross sections (S-Matrix elements) for particles produced at accelerators. 
In the early 80's the idea of the Inflationary universe was introduced by Guth, Linde and 
Starobinsky [1] and the mechanism for inflation was a scalar field (the inflaton) evolving 
in time in a background (time dependent) gravitational field. The initial calculations were 
done using a classical scalar field. It was natural for people to then try to understand how 
to solve time evolution problems in the quantum domain. In another area of Cosmological 
importance, Emil Mottola looked at the problem of particle production of "free" (apart from 
Gravity) scalar mesons in De Sitter Space as a possible mechanism for the solution of the 
cosmological constant problem. [2] . These two problems were the basis for renewed interest 
in the problem of how one solves numerically initial value and backreaction problems in 
quantum field theory. At that time even the correct formalism for doing this was not well 
understood by particle physicists. In trying to solve this problem, one has the choice of 
the Schrodinger picture (functional Schrodinger equation) or the Heisenberg picture (path 
integral or Schwinger- Dyson equations). Even though, in the Heisenberg picture, Schwinger 
and Keldysh [3] had introduced the correct formalism for treating this problem, the closed 
time path formalism (CTP), at that time very few relativistic field theory calculations had 
been attempted and the question of how to do practical calculations which required renor- 
malization had not been addressed. Also at that time it was not realized that if one used 
perturbation theory in the CTP formalism, that the resulting perturbation series (beyond 
leading order) was secular [I] . 

The first study of inflation which included quantum effects was In 1985, when Guth and Pi 
[5] studied the problem of the quantum roll from the top of an inverted harmonic oscillator in 
the free field limit to see to what extent classical ideas on inflation were modified by quantum 
effects. The formalism used to study this problem was the functional Schrodinger equation 
which was the generalization of the ordinary Schrodinger equation to quantum field theory. 
This study however did not require renormalization. The functional Schrodinger equation 
could be derived by the Dirac action principle. This then allowed people to study time the 
quantum evolution process using gaussian trial wave functions (i.e. variational approxima- 
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tions). This approximation to the field theory was related to solving the field theory in 
the Heisenberg picture in a mean field (Hartree approximation or leading order in large-N 
approximation). The first conceptual problems that needed to be solved in the scalar field 
theory in this approximation were how to disentangle the infiinities that arose from choosing 
(potentially unphysical) initial data from those that were truly effects of renormalization. 
These issues and their solution was clarified in the work of Cooper and Mottola and Sami- 
ullah and Pi [6J. The initial value problem was solved in the gaussian (leading order in 
large-N) approximation by choosing initial states that corresponded to finite energy density 
and number density with respect to the vacuum of the gaussian approximation. The renor- 
malization issues were first understood by doing a WKB analysis of an adiabatic expansion 
of the Green's functions. This approach was related to the method of adiabatic regulation 
used to study free fields in background gravitational fields [7] . Later it was shown that one 
could study the renormalization group flow of the coupling constant with momentum and 
verify that a more standard renormalization which looked for the discrete momentum space 
regulated answer converging to the continuum renormalization being the simplest approach 
[8]. Although the functional Schrodinger equation coupled with variational ansatzes led to 
approximations which conserved probability and energyunlike what happens in the case of 
variational approximations in quantum mechanics, it was difficult to find ways to go beyond 
the Gaussian approximation in a systematic fashion. On the other hand the CTP formalism 
has a nice path integral representation which allowed for ordinary perturbation expansions 
as well as 1/N expansions which were systematic. The downside of the CTP formalism was 
that the latter two expansions directly applied to the path integral turned out to be secu- 
lar and did not exactly conserve probabilities. This serious defect would later be cured by 
looking at expansions based on the generating functional of two particle irreducible graphs. 

It was in order to better understand back reaction as a solution of the cosmological 
constant problem that Cooper and Mottola first studied as a "toy" model, backreaction 
in the Electric Field case. [9]. The degradation of the electric field as a result of pair 
production was the test case for the degradation of the cosmological constant as a result of 
pair production. While this pair production problem was being investigated, there was an 
independent compelling reason for studying this problem. In order to understand particle 
production following relativistic heavy ion collisions people began investigating in the context 
of transport theory, a flux tube model based on Schwinger pair production, which converted 
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this toy problem into one with experimental consequences [10J. In light of the renewed 
interest in particle production from semi-classical gluonic fields we have recently undertaken 
a study of the quantum back reaction problem in SU(3) QCD in 3+1 dimensions in the hope 
of seeing what one can learn about jet production at RHIC and LHC and also the initial 
gluon condensate state. 

In this talk we would first like to review what is known analytically about pair production 
from strong electric and chromoelectric fields, both static and time dependent. Then we will 
concentrate on making the problem more realistic by solving simultaneously as an initial 
value problem the Dirac equation in an external electric field, and the update equation for 
the electric field as it degrades and oscillates due to expansion and particle production. We 
will first review previous results on pair production and back reaction in QED mostly in 1+1 
dimension, where the kinematics relevant for Heavy Ion collisions were used. Next we will 
discuss current work on solving the back reaction and particle production problem ins SU(3) 
QCD in 3+1 dimensions. Finally we will discuss how one would go about adding interactions 
so that the issue of thermalization of the fields during expansion can be addressed. 



II. PAIR PRODUCTION FROM A STRONG ELECTRIC AND CHROMOELEC- 
TRIC FIELDS 

In order to pop a pair of fermions (or bosons) out of the vacuum using strong Electric 
Fields one must supply an energy eEx in a Compton wave length x ~ h/mc which is at 
least twice the rest energy of the pair ( 2mc 2 ) . The critical value of the electric field for 
this to happen is of order 

eEh/mc = 2mc 2 (2.1) 

or 

Ecriucai ~ 2m 2 c 3 /eh. (2.2) 

This critical electric field is not yet attainable using lasers, but the analogue process of 
producing chromoelectric fields by colliding heavy ions does lead to initial energy densities 
that are above the critical value. The dimensionless variable relevant for this process is thus 

(2.3) 



Ecriucai 2ehm 2 c 3 
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In order to see how this variable arises, one can make the following simple tunneling picture 
[TT] of the non-perturbative process of pair production. One imagines that one has an 
electron bound in a potential well of order |Vq| « 2mc 2 and one then applies a constant 
electric field which leads to a one dimensional extra potential of eEx to the (say square 
well) potential of depth 2mc 2 . On then finds that the ionization probability is proportional 
to the WKB barrier penetration factor: 

^V /eE 



exp 



rv /ea 

-2 / dx (2m(Vb - \eE\x) 1/2 ) 
Jo 



Sm 2 c J /3ehE) 



(2.4) 



In what follows we will set h = l;c = 1. We see that when E > E cr i t i ca i there is no 
exponential suppression of pair production. A more careful calculation discussed below 
involves determining the imaginary part of the vacuum persistence function in the presence 
of the external field. 



A. Constant electric and chromoelectric field results 

In his classic paper in 1951 Schwinger derived the following one-loop non-perturbative 
formula 

dW _ e 2 ^ 1_ 

71=1 

for the probability of e + e~ pair production per unit time per unit volume from a constant 
electric field E via vacuum polarization [12J by using a proper time method. In case of 
charged scalar field theory the corresponding result is given by 

dW e 2 E 2 4^ (-l) n+1 ^ t 

e"W. (2.6) 



d A x 8vr 3 ^ n 2 

n=l 

The result of Schwinger was extended to QCD by Claudon, Yildiz and Cox [13] . However 
the px distribution of the e + (or e~) production, dW / (d A xd 2 pT), could not be obtained by 
using proper time method of Schwinger. A WKB approximate method was used for this 
purpose by Casher et. al. [H], but an exact method to do this problem (of determining 
the transverse distribution of pairs) was not found until recently [13]. For QED the WKB 
analysis gave the correct answer which depended only on the energy density of the electric 
field. However for QCD, the WKB answer was similar to the exact answer for QED but 
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incorrect for QCD in that the exact answer depended on both Casimir invariants of SU(3). 
In the case of fermions in QED one finds for the transverse distribution of fermion pairs: 

dW \eE\ 



L Log[l 



d 4 xd 2 pT 47T 3 

The corresponding result for the charged scalar production is given by 



(2.7) 



dW 



eE\ _ Pr+ vn 
-Log[l + e 71 ~T^T 



in 



(2.8) 



d A xd 2 pT 

In QCD the transverse distribution instead depends on two independent Casimir invari- 
ants of SU(3): Cx = [E a E a ] and C 2 = [d abc E a E b E c ] 2 , where E a is the constant chromo- 
electric field with color index a=l,2,..8 |15j . Nayak obtained the following formula for the 
number of non-perturbative quarks (antiquarks) produced per unit time, per unit volume 
and per unit transverse momentum from a given constant chromo-electric field E a 



dN, 



<?.<? 



dtd 3 xd 2 pq 



4tt 



1 3 7r(p|,+m 2 ) 



(2.9) 



where m is the mass of the quark. This result is gauge invariant because it depends on the 
following gauge invariant eigenvalues 



Ai 
A 2 
A 3 




- 
Cx 





cos9 , 

cos (2tt/3 - 6) 
cos (2tt/3 + 0) 



(2.10) 



where 9 is given by 



cos 2 3^ =3C 2 /Cf. 



(2.11) 



The integration over in eq. (2.9) reproduces Schwinger's proper time result, extended 



to QCD, for total production rate dN/d 4 x [13J. The exact result in eq. (2.9) can be 
contrasted with the following formula obtained by the WKB tunneling method [H] 



dN. 



q.q 



■\gE\ 



ln[l 



(2.12) 



dtd 3 xd 2 pr 47r 3 

which does not reproduce the correct result for the pr distribution of the quark (antiquark) 
production rate from a constant chromo-electric field E a . For soft gluon production Nayak 



and van Nieuwenhuizen [15j found in the Feynman T'Hooft gauge [T6] 

5^ - ^ £ 1^1 mp + rw]. (2.13) 

This was shown to be independent of the gauge fixing choice by Cooper and Nayak |17j . 
B. Schwinger pair production in SU(3) gauge theory 

In the background field method of QCD the gauge field is the sum of a classical back- 
ground field and the quantum gluon field: 

A% - A% + Ql (2.14) 

where in the right hand side A a is the classical background field and is the quantum 
gluon field. The gauge field Lagrangian density is given by 

£ g au gc = -\f%[A + Q] F a ^[A + Q]. (2.15) 

The background gauge fixing is given by by [16] 



D lt [A]Q> ta = 0, (2.16) 

where the covariant derivative is defined by 

Df[A] = 5 a % + gf abc A^. (2.17) 

The gauge fixing Lagrangian density is 

C gi = -^[D,[A]Q^] 2 (2.18) 

where a is any arbitrary gauge parameter, and the corresponding ghost contribution is given 
by 

£ gh ost = X a Df[A] D^ bc [A + Q] x c = x a K ab [A, Q] x b • (2.19) 
Now adding Eqs. (2.15), ( 2.18[ ), and (2.19), we get the Langrangian density for gluons 



interacting with a classical background field: 

£ g iuo„ = -\F« v [A + Q]F a ^[A + Q] - i- [D,[A}Q^] 2 -x a K ab [A,Q] X b . (2.20) 



To discuss gluon pair production at the one-loop level on considers just the part of this 
Lagrangian which is quadratic in quantum fields. This quadratic Lagrangian is invariant 
under a restricted class of gauge transformations. The quadratic Lagrangian for a pair of 
gluon interacting with background field A ^ is given by 

C g% = l -Q» a Mf v [A]Q vb (2-21) 

where 

M;l[A] = r] flu [Dp(A)D p (A)] ab - 2gf abc F* v + ( ± - 1 ) [D^A) D u (A)] ab (2.22) 

with 77^ = (-1, +1, +1, +1). 
For our purpose we write 

M£[A) = M^ a=l [A) + a' \D,{A) D M)] ah (2.23) 
where a' = (- — 1). The matrix elements for the gauge parameter a—1 is given by 

M$., a=1 [A] = VfMU [D p (A)D p (A)] ab - 2gf"* c F* v (2.24) 
which was studied in [15] . In this approximation, the ghost Lagrangian density is given by 

£ g host = X a D; b [A}D^[A] X c = X a K ab [A] X b (2.25) 

The vacuum-to-vacuum transition amplitude in pure gauge theory in the presence of a 
background field A® is given by: 

+ ( I ) A = J [dQ] [d X ] W] e^ s+s ^ +s ^ . (2.26) 

For the gluon pair part this can be written by 

, / n ) = = J = e lS °« 2 27) 

+ K 1 Z[0] J[rfg] e i/^Q^^[o]Q" 6 1 ' 

where is the one-loop effective action. The non-perturbative real gluon production is 
related to the imaginary part of the effective action which is physically due to the 
instability of the QCD vacuum in the presence of the background field. The above equation 
can be written as 



Z[0] Det-V2M$[0] 



< I ) A _ = ^ = ^_ U9 n Z J = (2.28) 
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which gives 

S °* = - iL ^!^m-^ = ~2 T ^K b M] - LnM-[0] ] . (2.29) 

The trace contains an integration over d A x and a sum over color and Lorentz indices. To 
the above action, we need to add the ghost action. The ghost action is gauge independent 
and eliminates the unphysical gluon degrees of freedom. The one-loop effective action for 
the ghost in the background field is given by 



s $Lt = ~iLn{Det K) = -iTr I — e ls lKlU|+ie| - e ls ^ A|+ * eJ (2.30) 

Jo s L 

where i^ ab [A] is given by (2.25). Since the total action is the sum of the gluon and ghost 



ds 



actions, the gauge parameter dependent part proportional to ( 1/a — 1 ) can be evaluated as 
an addition to the a = 1. 

The non-perturbative gluon pair production per unit volume per unit time is related to 
the imaginary part of this effective action via 

™ lmC ^^. (2.31) 



dtd 3 x d 4 x 
This is the general formulation of Schwinger mechanism in pure gauge theory where M®£ [A] 



is given by eq. (2.23) and AA] is given by eq. (2.24). In [15] this expression was 

evaluated for a = 1 and the final expression for the number of non-perturbative gluon (pair) 
production per unit time per unit volume and per unit transverse momentum from constant 
chromo-electric field E a is given by [To] 



dNg^g 1 ^ , - , . r. 



J2 \9^j\ Ml + e i^i]. (2.32) 



dtd 3 xd 2 pT 47r 3 | 

After this calculation was done in a = 1 gauge, Cooper and Nayak showed by explicit 
evaluation of the extra term proportion to a — 1 that the result for the particle production 
rate was independent of the gauge fixing parameter a [17j . 

Recently there has also been some progress to extending this result to time dependent 
fields. Using a formal operator shift theorem [18], Nayak and Cooper were able to show: 



dW \eE(t)\ T r ^T +m \ 

7l -Log[l + e * ™ ]. (2.33) 



d 4 xd 2 pT 8n 3 
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For Fermion pair production we get instead: 



dW \eE(t)\ T ri v\+ m \ 

- ' Log[l -e ]. (2.34) 



d 4 xd 2 p T 47r 3 

These results came from evaluating the one loop Action. For the Boson case we obtained 

S$ = —J C °-fd*x[ d 2 p T e^^[- _ eE(t) 
3 16tt 3 J s J J s sinh(se£(t)) J v ; 

wheras in the fermion case we obtained instead: 

S w = -At I" — I d A x ! d 2 p T e ls{p T +m2+ie \eE(t) coth(seE(t)) - -]. (2.36) 
8vr 3 J s J J s 

The way the above results were obtained was to start with the effective action for scalar 
field theory where 

d 

M[A] = (p-eA) 2 -m 2 ; ^ = 2 — 

The Action is 

= ^Trln[(p — eA) 2 — m 2 ] — iTr\n[p 2 — m 2 ] 

= -i J" ^ j d A X < x |[ e -»[(P-^) 2 -m 2 -«] _ e -i8(jP-m*-ie)^ x > 

Choosing Axial gauge A 3 = and A = —E(t)z 

roc r+oo r+oo r+oo r+oo 

S {1) = - i — dt{t\ dx(x\ dy(y\ dz 

Jo & J — OO J— oo J — oo J — oo 

^| ^ e -is[{p +eE{t)zY-pl-f? T -rn*-ie] _ ^(pW-ie ) } \ Z ) \y) \ X ) \t) (2.37) 

Inserting complete set of \pr) states (J d 2 pT\pr) {pt\ = 1) and using (q\p) = ^^e iqp we 
obtain 

dt(t\ dz(z\e- ls[ (- i ^ +eE{t)z)2 - p ^z)\t)- dt dz—j. (2.38) 

) */ — OO V V 



This expression contains the noncommuting quantities E(t) and d/dt. To evaluate these 
terms we derived a new shift theorem for operators: 



/ 

J — c 



dx(x\ e - [{a{y)x+h ^ )2+b ^ )+c{y)] \x) f(y) 

= [ + °°dx(x e-^* 2+b ^ +c ^ x -J-±)f( y ). (2.39) 

J-oo X a(y) dy a(y) dy I 
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This Shift Theorem also implies the interesting result. 

/ +0 ° dx e - {f{y)x+ i )2 g{y) = [ + °° dx e~^ x2 g(y) = ^ (2.40) 

i-oo J-oo f{y) 

The results we found for the time dependent field have the remarkable feature that they are 
equivalent to Schwinger's original expressions for the effective action with the substitution 
E — > E(t). That is the adiabatic approximation appears to give the exact result for the 
action. Although this is initially surprising, it is not without precedent. A related result for 
the one loop effective action was found recently by Fried and Woodard [2U] , using Fradkin's 
formulation of the path integral, for the case of an electric field pointing in the z direction 
which arbitrarily depends on the light cone time coordinate x + = (x° + a; 3 ). Explicitly they 
found that the action integrated over momentum was also equivalent to the adiabatic result 
in the variable x + , namely for the fermion action they obtained: 

T 1 [A] = -iL[A} = ^ J d4x J" ^e- lsm2 {eE(x + )s coth(eE(x + )s) -l} . (2.41) 

III. PRODUCTION AND TIME EVOLUTION OF A QUARK- ANTIQU ARK 
PLASMA 

Our model for the production of the quark-gluon plasma begins with the creation of a 
flux tube containing a strong color electric field. If the energy density of the chromoelectric 
field gets high enough as discussed earlier, the quark-anti quark pairs (or gluon pairs) can be 
popped out of the vacuum by tunneling. For simplicity, first we discuss pair production (such 
as electron-positron pairs) from an abelian Electric Field and the subsequent quantum back- 
reaction on the Electric Field. We will give extensive results for 1 + 1 dimensions [21] . In the 
next section we will discuss extending this work to QCD in 3+1 dimension. Our initial value 
problem will be that the Electric (or Chromelectric) field is given at an initial (proper) time, 
and the fermions initially are in the vacuum state (no particles initially present). Having 
a semi-classical Electric field interacting with fully quantum fermions is the first term in a 
large-N expansion of N-flavored QED. [22] . We assume as a reasonable first approximation, 
guided by hydrodynamical considerations, that the kinematics of ultrarelativistic high energy 
collisions results in a boost invariant dynamics [23] in the longitudinal (z) direction (here z 
corresponds to the axis of the initial collision) so that all expectation values are functions of 
the proper time r = yt 2 — z 2 .We introduce the light cone variables r and r/, which will be 
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identified later with fluid proper time and rapidity . These coordinates are defined in terms 
of the ordinary lab-frame Minkowski time t and coordinate along the beam direction z by 



z = t sinh 77 , t — T cosh r\ . 
The Minkowski line element in these coordinates has the form 

ds 2 = -dr 2 + dx 2 + dy 2 + r 2 dif . 
Hence the metric tensor is given by 

= diag(-l, 1, l,r 2 ). 
The QED action in curvilinear coordinates is: 

S = J d d+l x{AetV)[- l -^^V^ 



where 



V^= (d^ + F^-ieA^ 
Varying the action leads to the Heisenberg field equation: 

(7^V M + m) * = , 

7° (d T + ^j+7i'9i + y(^- ieA ri ) + m 
and the Maxwell equation: E = E z (t) = —A v (t) 



* = 0, 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



We expand the fermion field in terms of Fourier modes at fixed proper time: r, 

*(*) = f[dk]J2[bs^X(r)e^e^ 

J s 

+c4(-k)^: k >)e- i Sr ip ' x ]. (3.9) 

The ip^ s then obey 

_ / r\ 1 \ 1 

(3.10) 



H(r) = 0, 
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Squaring the Dirac equation: 



-7 ( — + — ) - «7_l • k ± - ry 7r T; + m 



X^- (3.H) 



iWxs = KXs (3.12) 
with \ s — 1 for s — 1, 2 and A s = — 1 for s = 3, 4, we then get the mode equation: 

(j^ + "k " /£to=0, (3.13) 

^ = < + ki + m 2 ; vr, = ^^. (3.14) 
The back-reaction equation in terms of the modes is 

^^ = -^E/Mk](ki + m 2 )A s | /l i| 2 , (3.15) 

s=l 

A typical proper time evolution of E and j is shown in fig. 4. Here an initial value of 
E = 4 was chosen, and here E is the dimensionless E/E critica i relevant for unsuppressed pair 
production. 

A. Spectrum of particles 

To determine the number of particles produced one needs to introduce the adiabatic bases 
for the fields: 

J s 

+4 0)t (-k; r)^ ks (r)e^ Wkc!r ]e ik - x . (3.16) 

The operators 6 s (k) and &£°\k; r) are related by a Bogolyubov transformation: 

&<°>(k;r) = ^aj r (r)6.(k)+^ r (r)dj(-k) 
4 0) (-k;r) = £/£(T)&.(k) + aE(T)dJ(-k) (3.17) 

One finds that the interpolating phase space number density for the number of particles (or 
antiparticles) present per unit phase space volume at time r is given by: 

n(k;r) = £ <0 in |&<°>t (k; r)6 (o) (k; r) | 0m ) = £ |^,( r )P (3.18) 

r=l,2 s,r 
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This is an adiabatic invariant of the Hamiltonian dynamics governing the time evolution 
of the one and two point functions, and is therefore the logical choice as the interpolating 
particle number operator. At r = r it is equal to our initial number operator. If at later 
times one reaches the out regime because of the decrease in energy density due to expansion 
it becomes the usual out state phase space number density. 

The phase space distribution of particles (or antipartcles) in light cone variables is 

d G N 

MT) = fik ^ T) = ^dxldkld V dk v (3 ' 19) 
We now need to relate this quantity to the spectra of electrons and positrons produced by 
the strong electric field. We introduce the particle rapidity y and m± = \fk\ + m 2 defined 
by the particle 4-momentum in the center of mass coordinate system 

k/i = (m± cosh y, k±, m± sinh y) (3.20) 

The boost that takes one from the center of mass coordinates to the comoving frame where 
the energy momentum tensor is diagonal is given by tanhr/ = v = z/t, so that one can define 
the "fluid" 4-velocity in the center of mass frame as 

= (cosh 77, 0, 0, sinh 77) (3.21) 

We then find that the variable 



k 2 

co fc = ymi + -| = A;% (3.22) 

has the meaning of the energy of the particle in the comoving frame. The momenta k v that 
enters into the adiabatic phase space number density is one of two momenta canonical to 
the variables defined by the coordinate transformation to light cone variables. Namely the 
variables 

v ; ' 2 \t-z 

have as their canonical momenta 

k T = Et/r - k z z/r k v = -Ez + tk z . (3.23) 

To show this we consider the metric ds 2 = dr 2 — T 2 drf and the free Lagrangian 

m dx 11 dx v ,„„,, 

i= 2 9 -lfad7 (3 - 24) 
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Then we obtain for example 

dr r/dr dt dr dz 

K = m Ts =m[ ^ )z ds + ^ )t d^ 

Et — k z z , „ , 
= — = A;% (3.25) 

The interpolating phase-space density / of particles depends on k v , kj_, r, and is 77- 
independent. In order to obtain the physical particle rapidity and transverse momentum 
distribution, we change variables from (77, k v ) to (z, y) at a fixed r where y is the particle 
rapidity. We have 

j Txdz dx]_ J f(k v ,k±,r) (3.26) 



d 3 N d 3 N 
E- 



d 3 k T[dydk\ 
where the Jacobian J is evaluated at a fixed proper time r 



J 



dk v /dy dk v /dz 
drj/dy drj/dz 



(3.27) 



m_|_ cosier? — y) dk v 



J= ' " V/ yy =^lr- (3-28) 
cosh 77 



We also have 



/c r = mj_ cosh(r] — y); k v = — rm± sinh(?] — y) . (3.29) 

Calling the integration over the transverse dimension the effective transverse size of the 
colliding ions A± we then obtain that: 

J^kJ = A± I ^' T) " J^kJ (3 - 30) 

This quantity is independent of y which is a consequence of the assumed boost invariance. 
Note that we have proven using the property of the Jacobean, that the distribution of 
particles in partical rapidity is the same as the distribution of particles in fluid rapidity 
verifying that in the boost-invariant regime that Landau's intuition was correct [21]. 

We now want to motivate the Cooper- Frye formula used to calculate particle spectrum 
in hydrodynamical models of particle production [25] . First we note that the interpolating 
number density depends on k v and kj_ only through the combination: 



u; k =\lml + ^ = k^ (3.31) 
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Thus f(kjj, k±) = f(k fl u fl ) and so it depends on exactly the same variable as the comoving 
thermal distribution! We also have that a constant r surface (which is the freeze out surface 
of Landau) is parametrized as: 

da 11 = A± (dz, 0, 0, dt) = A ± dr](cosh rj, 0, 0, sinh 77) (3.32) 

We therefore find 

k^da^ = A±m±r cosh(r/ — y) = A±\dk r) \ (3.33) 

Thus we can rewrite our expression for the field theory particle spectra as 

d?N r r 

ndydk2 =A ± j dk v f(k v , k ± , r) = J /(/*%, r)PAx M (3.34) 

where in the second integration we keep y and r fixed. Thus with the replacement of the 
thermal single particle distribution by the interpolating number operator, we get via the 
coordinate transformation to the center of mass frame the Cooper-Frye formula. 

The boost invariant assumption leads to an Energy Momentum tensor which is diagonal 
in the(r, rj, x±) coordinate system which is thus a comoving one. In that system one has: 

= diagonal {e(r), Pll (r),p ± (r),p ± (r)} (3.35) 

We thus find in this approximation that there are two separate pressures, one in the longitu- 
dinal direction and one in the transverse direction which is quite different from the thermal 
equilibrium case. However only the longitudinal pressure enters into the "entropy" equation, 

e + p\\=Ts (3.36) 
d(er) „ 

d(sr) _ Ej eta 
dr ~ T 

In the out regime we find as in the Landau Model that st = constant. The energy density 
as a function of proper time is shown in Fig. [TJ For a one-dimensional boost invariant flow 
we find that the energy in a bin of fluid rapidity is just: 



dE 
dr] 



J T^da^ = A ± rcoshr]6(T) (3.37) 



which is just the (1 + 1) dimensional hydrodynamical result. Here however e is obtained 
by solving the field theory equation rather than using an ultrarelativistic equation of state. 

17 



20 



10 




20 



40 



60 



80 



100 



120 



FIG. 1: Time evolution of et as a function of r. 

This result does not depend on any assumptions of thermalization. We can ask if we can 
directly calculate the particle rapidity distribution from the ansatz that we divide the energy 
in a bin of fluid rapidity by the energy of a particle assuming that the fluidity rapidity is 
equal to the particle rapidity? 

dN 1 dE A i 



drj m cosh 77 drj m ^ ^ ^ 

We see from Fig. [2] that this works well even in our case where we have ignored rescattering, 
so that one does not have an equilibrium equation of state. In the field theory calculation the 
expectation value of the stress tensor must be renormalized since the electric field undergoes 
charge renormalization. We can determine the two pressures and the energy density as a 
function of r. Explicitly we have in the fermion case. 



e t 



(Trr) 



rS s / [dk]Rrr(k) + E 2 R /2 , 



where 



R TT (k) = 2(p 2 ± + m 2 )(s+|/ + | 2 - g \f-\ 2 ) - u 
-(p 2 ± +m 2 )(7r + eA) 2 /(8uj 5 T 2 ) 



Pll (r)r 2 =< T m >= r 



S s J [dkjXsirRnnik) - l -E 2 R r 2 



(3.39) 
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FIG. 2: The ratio of the approximate rapidity distribution £'/(mcoshy) {dN/dy) and exact ra- 
pidity distribution as a function of r. 



where 



and 



R m {k) = 2\f + \ 2 - (2u)-\u + Asvr)- 1 - A s ei/8a/V 2 

-X s eE/8uj 5 - X s n/4u 5 T 2 + 5n\ s (n + eA) 2 / {IQu 7 t 2 ) 



P±{r) =< T yy >=< T xx > 

= ( 4r )" 1 E / WKpKpI + m'r'Rrr - 2\irplR vv } 

s J 

+E 2 R /2. 



(3.40) 



(3.41) 



Thus we are able to numerically determine the dynamical equation of state pi = Pi(e) as a 
function of r. A typical result is shown in Fig. [3j 



IV. QCD BACK-REACTION PROBLEM WITH CYLINDRICAL SYMMETRY 



Previously we have discussed the constant Chromoelectric field problem and have shown 
that the answer for gluon production is independent of the Gauge Fixing term in covariant 
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FIG. 3: Proper time evolution ofp/e. 

gauges. The answer depended on two Casimirs. If the Schwinger mechhanism is important 
for these experiments in producing the initial plasma, one might be able to find an experi- 
mentally observable effect at RHIC and LHC for those jets coming from the semi-classical 
gluon plasma produced at RHIC following the collision of heavy Ions. One possibility is that 
event by event the transverse distribution of jets produced following a heavy ion collision 
might depend on the values of both Casimirs and not just the initial energy density present 
in the semi-classical gluonic field. Here we will only consider pair production of quarks, 
described by a field ip(x) and satisfying Dirac's equation: 

[r(d^-gA^x))+m]ij(x) = 0, (4.1) 

interacting with a classical Yang-Mills field A^(x) = A^(x) T a , where T a are the generators 
of the 577(3) algebra, and satisfying a back-reaction equation given by: 

DfF b ^{x) = g ( [4>(x),r(x)T a 4>(x)} >/2 , (4.2) 

with Df = 5 ab d„ + g f abc Afa). 

We want to consider cases when the plasma starts out in a hot equilibrium state with 
cylindrical symmetry, and then expands into the vacuum. We introduce the fluid rapidity 
(77) and proper time r coordinates which are components of aT 4 = (r,p,#, 77), and which are 
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related to Cartesian coordinates by: 



t = t cosh f] , z = t sinh 77 , x = p cos # , U = P sin # . 



(4.3) 



For a boost invariant expansion, the classical gauge fields are restricted to be in the re- 
direction and depend only on r. We also consider only the a = 3 and a = 8 gauge fields 
which carry all colors. This will allow us to explore the entire Casimir space. Then, using 
the Gell-Mann representation for the A" matrices, 



/ 



\ 



\ 



A 3 v (r) + Al(r)/V3 

-^(t)+A;(t)/V3 

-2A 8 v (t)/V3J 

The two Casimir invariants for £77(3) are given by: 

b T?C 1 2 



(4.4) 



Ci = £ a £ a , and C 2 = [ d ahc E a E b E c ] 



(4.5) 



where d abc are the symmetric SU (3) structure factors. Notice that F* F a '^ u is a Lorentz 



invariant and equal to C\. In our coordinate system, Eq. (4.2) becomes: 

ldE a g 



t Or 



(4.6) 



Eqs. (4.1) and (4.6) are the equations we want to solve. 



For the fermi field operator, we put: 



^(x) = S(6 1 r ] )^(x)/V^, 
where S(9,rj) = S p (8) S T (r]) is a Lorentz transformation given by: 



(4.7) 



S p (0) = exp[0 7 y/2 



S t (t]) =exp[77 7 u 7 72 



(4.8) 



and where Six) now satisfies: 



[7^(r,p)V M -m]0(x)/v^ = O. 



(4.9) 



Here 7 m (t, p) are given by: 



7 T = 7°, fW, 7 e (p)= 7 2 /P, 7 ^(r)= 7 3 /r, 



(4.10) 
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and V M = <9 M + II M + A M (r) is the covariant derivative in the boost-invariant system: 



n e = 7V/2 , n„ = 7 V/2 • 



(4.11) 



In cylindrical coordinates the canonical fields obey: 

{ Mr, p, 6, rj), &( T| pf, 6', rj') } = 8 a>a , 5{P ^ W ~ *0 " V) 



'PP 



(4.12) 



which is just the usual anticommutation relation. We can now write the Dirac field operator 
in this curvilinear coordinate system by the expansion: 



00 dkr, f°° ki dk\ 



00 2tt j 



2vr 



E E 



/i=±l m=— 00 



x { <JL»(r, P, f , 17) + 4"!t, m ^ ^(r, P, M) } . (4-13) 



where: 



HO 



(4.14) 



with X — hk \ , and where /i = ±1. 



(fc) / /i\ 1 6 M 

V2 Ue ! ( m+1 ) fl 4 +1 (M 



(4.15) 



with eigenvalues A = and the helicity /i = ±1. Orthogonality is given by the relation: 



pdp / (WxC 1 L(p^)xl?i Ilf (p,e) = V'^,-'(27r) 



S(k ± - k' ± ) 
y/klk^ 



The mode functions 



^(±W( r ) satisfy: 
i<9 T + m -n kv (r) - ihk± 



( r ) + i<9 T - m / \ 4>f\ kri , k± I 



where we have defined 7Tfc (r) by: 



(4.16) 



(4.17) 



^(t) = (A^-(/A(r))/r. 



(4.18) 



Eq. (4.17) is the equation we want to solve numerically as a function of r for some given 



initial spinor at r = 1. 
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Now at r = 1, we will set A(l) = 0. This means that near t = 1, Eq. (4.17) becomes: 



id T + m —k v — ihkj 
—k n + ihk± id T — m 

which have positive and negative frequency solutions of the form: 



J 0(+);k v ,k ± [ 
; 0(-);fe,,fc_L 



(4.19) 



6 {h ' +) (t) 



'^>0;k v ,k ± — rn 



( 



2lu< 



0;k v ,kj 



\ 



T 



' UJ 0;k v ,k ± — m 
^ UJ 0;k v ,k ± 



/ k v — ihk± \ 

1 



V 



exp[-^ ;fc t) ,A ;± (T - 1)] 



expf-Mc^^r - 1)] 



(4.20a) 



(4.20b) 



where cj 0; ^ 



k± 



k 2 + k ± + m 2 . These solutions are orthogonal: 



a=± 



and complete: 



(h,A)* 



(r) = 5 Q 



■ 



A=±l 



So at t = 1, we choose our solutions of Eq. (4.17) so that: 

^{a);k n ,kS^ = ^0 (a);k v ,k ± U) ' 



(4.21) 



(4.22) 



(4.23) 



for a = ± and where A = ±1 labels the initial positive and negative frequency solutions of 



Eq. (4.19). The r-dependent solutions will then be numerically stepped out from the values 



at r = 1. 

Maxwell's equation becomes: 

d T E(r) = - 



dk^ 

,27 



k± dfcj 
2tt 



3k v ]k ± ( T ) 



(4.24) 



h=±l 



where ^ x (t) is given by the positive energy solutions of the Dirac equation only: 



<(>*>+) t/ 



(4.25) 



'fc,,*x in 



Here, 0^'t x ( r ) is the two-component positive energy spinor: 

dh,+) 



(-);kr,,k± 



J) 



(4.26) 
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and a x the Pauli matrix. Dirac's Eq. (4.17) and Maxwell's Eq. (4.24), are the update 



equations we want to solve simultaneously. 

The definition of interpolating particle number for the time evolution problem is some- 
what arbitrary. The main need is for it to interpolate between initial an final definitions of 
the number operator and for it to have some adiabatic properties. We choose to define them 
here in terms of the exact solutions of the Dirac equation in the absence of external fields. 
These zeroth order spinors are given by: 



dh,X) 



y O(+);fc 77 ,fc i 
,(M) 



(4.27) 



where 0otvfc±( r ) gi yen by Eqs. (4.20). These spinors are also orthogonal and complete. 
Expansion of the field operator in the zeroth order spinors then requires that the creation 
and annihilation operators ^.cwk m( r ) become time- dependent. That is: 



oo 2tt jo 



k± dfcj 
2tt 



+oo 



E E E ^ ± , m (r)^ ± , m (r,pAv), (4-28) 



h=±l A=±l m=— oo 



Because of the orthogonality of the initial spinors, we see that the AoVxmW operators 
obey the same commutation relations as the time- independent ones at equal time: 



i(h',A')t 



f ,S(k x -k' ± 



(4.29) 



and so we can use them to define number operators at time r. Again, it is traditional to 
define particle and anti-particle operators at time r by: 



J ^0;k v ,k ± ,m\ l ) 



b {h) (t) 



and 



A {h '- ] (t) 



ol { - h)] (r) 



(4.30) 



As before we can determine the adiabatic number operator from the Bogoliubov transforma- 
tion. The overlap between the adiabatic wave functions and the exact ones is : C^'^ \r) 
is given by: 

^V) = ^>)^>), ( 4 - 31 ) 

and is independent of m. So the creation and annihilation operators are related by the 
expression: 



i(M) (_ \ _ r (h-X\') ( , Uh,\>) 

J ^0;k v ,k ± , m \' ) ~ 2-~i k v ,k ± ^ > k v ,k ± ,rn ' 



A=± 
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which is a Bogoliubov transformation of the operators. Now let n^^ k (r) be the number 
of particles produced at time r per unit volume in a given mode and helicity h: 



d 6 iV( 



fc„fcx,^v ) d ^ fci dfc± dfj pdpde - 
n^ k± ^(r) is given by the formula: 



(4.33) 



(4.34) 



(ft) 



' (t) 

v ,k±,<p\ i 



d 6 A^(-i 



dfc^ A;j_ dfcj_ d0 d?7 p dp d# 



C 



(h;+,+), 



T 



(4.35) 



and is independent of 0. Explicitly, | C^'^ (t) | 2 is given by 



C 



(h;+,+). 



u k v ,k A 



rn 



2ujk 



r + 



k v + ih k± 

Ukr,,k x - m 



(-);k v ,kj 



(4.36) 



The current vanishes at r = 1, as required. 



V. FLAVOR LARGE-N EXPANSION IN A 2PI APPROACH 



It is important to know whether rescatterings due to interactions betweeen the quarks 
and gluons will change the final transverse distribution functions from those found in the 
leading order in flavor large-N. To include interactions among the quarks and gluons one 
would solve the coupled Schwinger-Dyson equations using the CTP formalism and a 2PI 
[25] effective action expanded in 1/N [27] ■ Here one would need to keep the background 
field formalism also to handle the background chromoelectric field. This formalism has been 
used in QCD to determine transport coefficients by Aarts and Resco [28J. The action for Nf 
identical fermion fields ip a (a = 1, . . . , Nf) then reads (We use here g^ u = diag(+, — , — , — )) 



S 



with 



and we use the notation 



B = >fD 



dx° / d 3 x 



x JC 
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+ S'gf + S'gh, 



ie 



(5.1) 
(5.2) 
(5.3) 



FIG. 4: NLO contribution to the 2PI effective action in the 1/Nf expansion. 

where C refers to the CTP contour in the complex-time plane. Note that we have rescaled 
the coupling constant with WNf, so that in the large Nf limit Nf goes to infinity while e 
remains finite after renormalization. To fix the gauge one can use a general linear gauge 
fixing condition, 

S g f = -^[/^] 2 - (5-4) 
Below we specialize to the generalized coulomb gauge: fo — 0, fi — <9j. The ghost part is 
not needed explicitly. 

The 2PI effective action is an effective action for the contour-ordered two-point functions 
D^fay) = (T c [A^x) A v {y)] ) , S ab (x,y) = (T c [M^Mv)] > ( 5 - 5 ) 
and can be written as 

T[S, D] = % - Tr[ ln[ D" 1 ] ] + \ Tr[ D 1 (D - D ) } 

- i Tr[ ln[ S- 1 ] ] - % Tr[ S \S -S )} + T 2 [S, D\ + ghosts , (5.6) 

where Dq 1 and Sq 1 are the free inverse propagators and contains all two particle irre- 
ducible vacuum graphs such as those in Fig. [4} 

For the gauge theory the next-to-leading order Schwinger-Dyson equations that result 
from varying the 2PI action are: 

S- 1 = Sq 1 - E , D- 1 = D l - II , 

with S the fermion and D the gauge field propagator. The self energies, depending on full 
propagators, are shown in Fig. |5j 

The back reaction equation is given by: 

V M F"" = (f) = -ig 2 Tr^S (5.7) 

To determine the interpolating number densities of quarks and antiquarks one can follow 
the procedure of Berges, Borsanyi and Serreau [29] and define these from the current. 
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FIG. 5: Self energies at next-to-leading order in large Nf QCD. 

Namely the associated 4-current for each given flavor is ~ ■ip'j^if). Fourier transforming 
with respect to spatial momenta, the expectation value of the latter can be written as 
Jj(t,p) = Tr[^S < (t, t,p) ], In terms of the equal-time two point function, its temporal and 
spatial components are 

J° f (t,p) = 2[l-2F*(t,t;p)], 
3 f (t,p) = -AF v {t,t;p). 

To obtain an effective particle number, Berges et. al. identify these expressions with the 
corresponding ones in a quasi-particle description with free-field expressions. These are given 
by 

J°/ qp \t,p) = 2[l + Q f (t,p)], 
ji QP \t,p) = -2[l-2N f (t,p)], 

where Qf(t,p) = rif — % is the difference between particle and antiparticle effective num- 
ber densities and Nj(t,p) = (n/ + %)/2 is their half-sum. The physical content of these 
expressions is simple: the temporal component J° directly represents the net charge density 
per mode Qf(t,p), whereas the spatial part J is the net current density per mode and is 
therefore sensitive to the sum of particle and antiparticle number densities. Identifying the 
above expressions, they define 

^Q f (t,p) = -F$(t,t;p), 
^-N f (t,p) = F v (t,t;p). 

Using these definitions and solving the backreaction problem to NLO in 2PI-1/N we would 
also be able to discover if there is time for the produced quarks and antiquarks to thermalize 
before hadronization time scale and to see if the constant field result for the transverse 
distribution will be modified by the interactions. This will be an important calculation to 
do in the future. 
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